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^ ■ Abstract. Equilibrium configurations of cold neutron stars near the minimum mass are studied, using the recent 

equation of state SLy, which describes in a unified, physically consistent manner, both the solid crust and the 
liquid core of neutron stars. Results are compared with those obtained using an older FPS equation of state of cold 
catalyzed matter. The value of M m i n ~ 0.09 Mq depends very weakly on the equation of state of cold catalyzed 
matter: it is 0.094 Mq for the SLy model, and 0.088 Mq for the FPS one. Central density at M m i n is significantly 
lower than the normal nuclear density: for the SLy equation of state we get central density 1.7 10 14 g cm -3 , to be 
compared with 2.3 10 14 g cm -3 obtained for the FPS one. Even at Mmin, neutron stars have a small liquid core 
of radius of about 4 km, containing some 2-3% of the stellar mass. Neutron stars with 0.09 Mq < M < 0.17 Mq 
are bound with respect to dispersed configuration of the hydrogen gas, but are unbound with respect to dispersed 
Oh 56 Fe. The effect of uniform rotation on the minimum-mass configuration of cold neutron stars is studied. Rotation 

increases the value of M m i n ; at rotation period of 10 ms the minimum mass of neutron stars increases to 0.13 Mq, 
5_l ■ and corresponds to the mass-shedding (Keplerian) configuration. In the case of the shortest observed rotation 

£jj | period of radio pulsars 1.56 ms, minimum mass of uniformly rotating cold neutron stars corresponds to the 

mass-shedding limit, and is found at 0.61 M for the SLy EOS and 0.54 M for the FPS EOS. 
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1. Introduction gest, that after crossing the M m ; n value, a neutron star 

will undergo an explosion (Blinnikov et al. 1984, Colpi et 

Mass of neutron stars is bound on the lower-density side al 1989; iggi) Sumiyoshi et aL 1998 ). 

by the minimum mass, M m j n . Equilibrium neutron star 

configurations with M < M min do not exist. Some astro- 0f course ' station and thermal effects can affect (ac- 
physical scenarios, in which a neutron star looses mat- tuall y " urease) the actual value of the minimum neu- 
ter, and reaches the critical value M min , were proposed tron star mass ' As we wl11 see ' station will also change 
(Blinnikov et al. 1984, Colpi et al. 1989, Colpi et al. 1991, the nature of the mmimum-mass configuration. In the 
Sumiyoshi et al. 1998). Basically, one considers a neutron P resent P a P er ' we wlU restnct ourselves to neutron stars 
star in a close binary system with a more compact, and bullt of cold cat alyzed matter. For newly born proto- 
more massive object (more massive neutron star, stellar neutron stars < both thermal and neutrino-trapping effects 
mass black hole). At a sufficiently small separation be- are lar S e > and are ex P ected to lncrese minimum mass to 
tween a neutron star and its companion, less massive, and °- 9 " L1 M ° (Goussard et al. 1998, Strobel et al. 1999). 
therefore larger, neutron star starts to lose mass due to However, M min for static, cold neutron stars still remains 
gravitational pull of its companion, the lost matter being the absolute lower bound on the mass of neutron star, 
accreted by the companion compact object. This process The calculation of M m ; n for neutron stars has a long 
is self-accelerating, because decrease of mass leads to in- history. The first correct estimate of M m ; n was obtained 
crease of neutron star size, making it even more suscepti- by Oppenheimer & Serber (1938). They showed, that 
ble to the mass loss. After crossing M min , no equilibrium the previous estimate obtained by Landau (1938) was 
configuration can be reached. Numerical simulations sug- much too small, and moreover based on incorrect argu- 

ments. Oppenheimer & Serber (1938) obtained an esti- 

Send offprint requests to: P. Haensel mate 0.17 M Q without taking into account nuclear interac- 
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tions. Then they argued, that inclusion of the contribution 
from the nucleon-nucleon interactions leads to a decrease 
of M m ; n to 0.03 — 0.10 M© , depending on the assumptions 
they made about the (very poorly known at that time) 
nucleon-nucleon interaction. It should be stressed, that 
their criterion of stability was based on energy arguments: 
a neutron star had to be stable with respect to transfor- 
mation into medium mass nuclei (they took calcium as an 
example). Calculation of Harrison et al. (1964) was based 
on the Harrison- Wheeler equation of state (Harrison et 
al. 1964), which described both the crust and the liquid 
core of neutron stars. From the present-day point of view, 
their equation of state above neutron drip was unrealistic. 
They obtained M m ; n = 0.18 M©, at a central density of 
only 3 10 13 g cm~ 3 . Using two versions of the Levingcr 
& Simmons (1961) baryon-baryon potential, Tsuruta & 
Cameron (1966) obtained M min = 0.11 - 0.13 M Q . Cohen 
& Cameron (1971), using the equation of state (EOS) 
based on the up-dated version of the Levinger-Simmons 
nucleon-nucleon potential (Langer et al. 1969), obtained 
-^min = 0.065 Mq. Finally, in their classical paper on the 
EOS and structure of neutron stars, Baym et al. (1971b, 
hereafter BPS), who used their own (BPS) EOS of the 
outer crust, and that of Baym et al. (1971a, hereafter 
BBP) for the inner crust and for the outer layers of the 
liquid core, obtained M min = 0.0925 M Q . 

The properties of the equilibrium configurations near 
M m ; n are sensitive to the EOS of neutron star mat- 
ter at subnuclear densities, i.e. for mass density signifi- 
cantly below normal (saturation) density of nuclear mat- 
ter, po — 2.7 10 14 g cm~ 3 (baryon density significantly 
below n = 0.16 fm -3 ). As we will show, a particularly 
important role is played by the EOS near the inner crust 
- liquid core boundary. Therefore, the EOS should give 
physically correct description of this transition region; a 
brutal ad hoc matching of the inner crust and liquid core 
are not sound for this purpose. The microscopic model, 
underlying the EOS, should be the same on both sides of 
the crust-core boundary, and should therefore be based on 
the same effective nuclear hamiltonian. 

In the present paper, we calculate M m ; n and study 
static configurations near minimum mass. We also study 
the effect of rotation on these configurations. We use 
two EOS, based on different effective nuclear hamilto- 
nians. These effective hamiltonians are: a recent SLy 
model (Chabanat et al. 1998, for applications see Douchin 
& Haensel 2000, Douchin et al. 2000, EOS in Douchin 
& Haensel 2001), and a somewhat older FPS model 
(Ravenhall & Pandharipande 1989, for applications see 
Lorenz et al. 1993). 

The plan of the paper is as follows. Equation of state 
near crust-liquid core interface is discussed in Sect. 2. In 
Sect. 3, we study equilibrium configurations near M m i n , 
and their dependence on the assumed EOS model. The 
problem of binding energy and stability of low-mass neu- 
tron stars is studied in some detail in Sect. 4. Possible ef- 
fects of elastic strain within the equilibrium configurations 
near M m - ln are estimated in Sect. 5. Effects of rotation are 
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Fig. 1. Comparison of the SLy and FPS EOS near the 
crust-liquid core transition. Thick solid line: inner crust 
with spherical nuclei. Dashed line corresponds to non- 
spherical ("exotic") nuclear shapes. Thin solid line: uni- 
form npe matter. Vertical dashes separate different dense 
matter phases. 

studied in Sect. 6. Finally, Sect. 7 contains discussion of 
our results and conclusions. 

2. Equation of state of cold catalyzed matter near 
the crust-core interface 

The properties of neutron stars close to M m ; n are deter- 
mined by the EOS of dense matter at subnuclear densities, 
and in particular in the density interval 0.1 po < p < po- 
At such densities, the actual structure of the ground state 
of matter results from the interplay of the surface and 
Coulomb terms in the total energy density. These tiny 
"finite size" terms, which are absent in uniform npe mat- 
ter, are difficult to calculate with high precision at the 
density of interest. However, these terms, very small com- 
pared to the main bulk energy contribution, decide what 
is the actual shape of nuclei, and whether or not matter 
with nuclear structures is energetically preferred over the 
spatially uniform npe matter. 

As we already stressed in Sect.l, it is mandatory to use 
the same effective nuclear hamiltonian to describe both 
nuclear structures (nuclei) and neutron gas in the inner 
crust, and the uniform npe matter of the liquid core. In 
the present paper we use a recently developed SLy model 
of neutron star matter (Douchin & Haensel 2000, Douchin 
et al. 2000, Douchin & Haensel 2001), based on the SLy 
effective nuclear interaction (Chabanat et al. 1997, 1998). 
For this model, nuclei in the ground state of neutron star 
matter remain spherical down to the bottom edge of the 
inner crust (Douchin & Haensel 2000). Transition to the 
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uniform npe plasma takes place at p e dge = 1.3 10 14 g cm~ 3 
(baryon density n c dgc = 0.078 fm~ 3 ). The crust-core phase 
transition is very weakly first-order, with a relative density 
jump of about 1%. 

Unfortunately, absence or presence of non-spherical 
(called also "exotic" or "unusual") nuclear shapes in the 
bottom layers of the inner crust depends on the assumed 
model of effective nuclear hamiltonian. As an example 
of a different effective nuclear hamiltonian, we will con- 
sider the FPS model (Pandharipande & Ravenhall 1989, 
Lorenz et al. 1993), which also gives a unified descrip- 
tion of both the inner crust and the liquid core. For the 
FPS model, the crust-liquid core transition takes place at 
higher density p c dgc = 1.6 10 14 g cm -3 (baryon density 
^edgc = 0.096 fm~ 3 ), and is preceded by a sequence of 
phase transitions between various nuclear shapes (Lorenz 
et al. 1993). The sequence of phase transitions starts at 
1.1 10 14 g cm -3 ~ ipo (baryon density rih = 0.064 fm -3 ), 
at which spherical nuclei immersed in neutron gas are 
replaced by the cylindrical ones. Then follows a transi- 
tion from cylindrical nuclei immersed in neutron gas to 
the slabs of nuclear matter interspaced with neutron gas, 
which at still higher density are replaced by nuclear mat- 
ter with cylindrical holes filled with neutron gas, replaced 
eventually by nuclear matter with spherical bubbles filled 
with neutron gas. This last phase persists up to p e dge- All 
phase transitions are very weakly first-order, with relative 
density jumps smaller than 1%. 

Both EOS are compared, in the relevant density in- 
terval, in Figjl|. In the vicinity of the crust-core interface, 
the SLy EOS is stiffer than the FPS one. Moreover, in the 
case of the SLy EOS the discontinuous stiffening (jump of 
compression modulus) at the crust-core transition is more 
pronounced than in the FPS case, where it has been pre- 
ceded by a sequence of phase transitions connected with 
changes of nuclear shapes. Such a sequence of phase tran- 
sitions implied therefore a gradual stiffening of matter. 

3. Mass, radius and central density of 
configuration close to M m ; n 

The mass-radius plot for the SLy EOS, in the vicinity of 
M m i n , is shown in Fig.^. The mass minimum is very flat, 
particularly on the lower-density branch. In contrast, the 
minimum of the mass-central density curve is very sharp, 
and the curve M(p c ) on the lower-density side is extremely 
steep, as seen in Figjj^. Configurations to the left of the 
minimum in Figj^ have dM/dp c < 0, and are therefore 
unstable with respect to small radial perturbations; they 
are not expected to exist in the Universe (see, e.g., Shapiro 
& Teukolsky 1983). 

Let us denote central density of the M m i n configuration 
by Pc,min- Configurations with p c ~ p c ,min are very loosely 
bound by gravitational forces (see below). Their radii are 
~ hundreds of kilometers, and a small difference in their 
mass is accompanied by a large difference in radii. For 
example, on a stable branch on Fig.^, a decrease in mass 
by 0.001 M Q (i.e., by one percent) implies an increase in 
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Fig. 2. Gravitational mass versus radius near M m ; n , for 
cold, static neutron stars. Calculation performed for the 
SLy EOS. Solid line - stable configurations, dotted lines - 
configurations unstable with respect to small radial per- 
turbations. Minimum mass configuration is indicated with 
a filled circle. 
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Fig. 3. Gravitational mass versus central density, in 
the vicinity of M m i n , for cold, static neutron stars. 
Calculations performed for the SLy and the FPS EOS. 
Solid line - stable configurations, dotted lines - configura- 
tions unstable with respect to small radial perturbations. 
Minimum mass configuration is indicated with a filled cir- 
cle. 



radius by ~ hundred kilometers (i.e., by more than a factor 
of two !). Simultaneously, as seen in Fig.|, for p c < pc. 

min 

(unstable branch) the increase of mass with decreasing 
central density becomes extremely steep. 

This behavior is related to the specific structure of neu- 
tron stars near M m ; n , and to the properties of the EOS of 
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Table 1. Configurations of minimum mass for cold, static neutron stars 

EOS M~ J c R M cote /M R^ e E^ d E^J~ d 
[M Q ] [10 14 g/cm 3 ] [km] [km] [10 51 erg] [10 51 erg] 

SLy 0.094 1.6 270 0.02 3.8 -1.0 0.5 

FPS 0.088 2.2 220 0.03 4.2 -1.2 0.2 



neutron star matter at subnuclear density. Neutron stars 
close to M m i n are loosely bound, huge spheres of the solid 
crust, of the radius of hundreds of kilometers, containing 
a tiny liquid core, which at M m i n constitutes 2% of the 
stellar mass for the SLy EOS, and 3% of mass in the case 
of the FPS model. 

In the context of the stability of equilibrium configura- 
tions, a particular role is played by the adiabatic index of 
the neutron star interior, T = (n\y/ P)(dP/&n\ } ). As we see 
in Figj^ the most dramatic steepening of the M(p c ) curve 
takes place at p c ~ 1.6 10 14 g cm' 3 for the SLy EOS, to be 
compared with p c ~ 2.2 10 14 g cm" 3 for the FPS EOS. At 
such central density, neutron stars have a tiny, stiff, liquid 
core with r core ~ 2. The average value of the adiabatic 
index within the crust is much lower, r cmst ~ 1.3, below 
the Newtonian threshold for instability with respect to 
small radial perturbation. Gravitational pull barely binds 
the equilibrium configurations, which are unstable with 
respect to small radial perturbations. 

At M m i n , the radius of the liquid core is 3.8 km for 
the SLy EOS, and 4.2 km for the FPS model. As the SLy 
EOS is stiffer for 10 13 g cm -3 < p < p than the FPS 
one, the SLy neutron stars at -M m ; n are less compact, and 
contain smaller and less massive liquid core. Parameters 
of neutron star models with minimum mass are collected 
in Table 1. 

For p c < /d c ,min, both the mass of the liquid core and its 
radius decrease very slowly with decreasing p c , while the 
total mass and total radius of the star increases extremely 
steeply. 

4. Binding energy of low-mass neutron stars 

Binding energy of a neutron star, -Ebind, is defined as the 
mass defect with respect to dispersed configuration of the 
same number of baryons, multiplied by c 2 . Dispersed con- 
figuration is characterized by negligible pressure and grav- 
itational interactions. Equivalently, one may define -Ebmd 
as a net work needed to transform a neutron star into a 
dispersed configuration. 

One may contemplate several definitions of dispersed 
configuration. Let us first consider the net work needed to 
transform neutron star, consisting of A baryons, into a dis- 
persed configuration, under the condition of keeping mat- 
ter always in the ground state. In this way, dispersed con- 
figuration will be a pressureless cloud of the dust of 56 Fe, 
with mass per nucleon mp c = mass of 56 Fe atom/56 = 
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Fig. 4. Gravitational mass versus central density, and 
binding energy (relative to dispersed 56 Fe - labeled by Fe, 
and relative to hydrogen gas - labeled by H) versus gravi- 
tational mass, for cold, static neutron stars. Calculations 
performed for the SLy and the FPS EOS. Solid line - sta- 
ble configurations, dotted lines - configurations unstable 
with respect to small radial perturbations. Minimum mass 
configuration is indicated with a filled circle. 

1.6587 10~ 24 g. Binding energy with respect to such a 
dispersed configuration is 

4fn C d=(^Fc~M)c 2 . (1) 

Another dispersed configuration is that of a pressureless 
hydrogen cloud. Binding energy in that case is 

^fn } d=(^H-M) C 2 , (2) 

where the mass of the hydrogen atom mjj = 
1.6735 10" 24 g. 

Binding energy of low-mass, cold, static neutron stars, 
versus stellar mass, in the vicinity of M m - m , is plotted in 
Fig. |]. As we see, neutron stars with M < 0.17 M Q are 
unbound with respect to 56 Fe. However, neutron stars are 
always bound with respect to the hydrogen cloud config- 
uration. At the same mass (or baryon number), unstable 



P. Haensel et al.: Minimum mass of cold neutron stars 



5 



configuration (dotted line in Figjj) is less bound than the 
stable one. 

The case of .Ebind < deserves a comment. Negative 
binding energy with respect to dispersed 56 Fe indicates, 
that neutron stars with M < 0.17 M Q are actually 
metastable with respect to transformation into a cloud of 
56 Fe dust. Their decompression (while locally keeping the 
matter in the ground state) would eventually lower the 
total energy, after overcoming an initial energy barrier 
resulting from gravitational binding of the outer layers. 
While being metastable with respect to large amplitude 
decompression, configurations with p c > p c ,min are, how- 
ever, stable with respect to small amplitude radial per- 
turbations; this is to be contrasted with an instability of 
configurations with p c < p c ,min- One expects, that equi- 
librium configurations with p c < p c .min, after being per- 
turbed, will eventually expand to huge white-dwarf like 
configurations with radii exceeding 10 4 km. 

5. Effect of elastic strains 

Configurations of hydrostatic equilibrium were calculated 
using the Oppenheimer-Volkoff equation (OV) of general 
relativity (see, e.g., Shapiro & Teukolsky 1983). The OV 
equation is derived under the assumption that matter is 
described by the ideal fluid stress tensor. However, at 
M ~ 0.1 M a neutron star is nearly all solid, except for a 
small liquid core containing only 2 — 3% of stellar mass. In 
contrast to liquid, a solid can sustain a shear strain, which 
contributes to the stress tensor. In the isotropic approxi- 
mation, the elastic shear-strain term in the stress tensor 
of the Coulomb crystal of the solid crust is determined 
by the shear modulus fi. The effect of elastic strain on 
the stellar structure can be roughly estimated using the 
model of Baym & Pines (1971). In the presence of elas- 
tic shear strain, static neutron star configuration can be 
non-spherical. Assuming for simplicity that the deviation 
from sphericity is quadrupolar, let us estimate possible 
cllipticity 8R/R, where SR is the absolute value of the 
difference between the polar and the equatorial radii. For 
a Coulomb crystal the maximum (i.e., breaking) strain 
does not exceed ~ 10~ 2 (Ruderman 1992). The Baym- 
Pines (1971) model predicts that 5R/R < 0.016 where 
the "rigidity parameter" b = (3 f v t ^dy/|-E grav |. Here, 
-Egrav is the gravitational energy of a neutron star and /3 
is a numerical factor ~ 1. For a 1.4 M neutron star, one 
typically has b ~ 10~ 5 , but for M ~ 0.1M Q one obtains 
b ~ 10" 2 (Carlini & Treves 1989). Therefore, we expect 
that the effect of elastic shear strain on the stellar radius 
at M ~ 0.1 M Q is 8R/R < 10~ 3 . 

6. Effects of rotation 

The shortest observed pulsar period is = 1.558 ms, 
which corresponds to rotation frequency = 641 Hz. 
For massive neutron stars of M <; 1.3 M Q , rotation at 
such frequency has a rather small offect on neutron star 
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Fig. 5. Gravitational mass, M, versus equatorial radius 
for low-mass neutron stars and effects of uniform rotation. 
For non-rotating neutron stars equatorial and polar radii 
are equal, i? oq = R. Solid line corresponds to non-rotating 
neutron-star models based on the SLy EOS. Short-dash 
lines correspond to uniformly rotating neutron-star mod- 
els based on the SLy EOS, the upper one for rotation fre- 
quency of the most rapid observed millisecond pulsar, and 
the lower one to rotation period of 10 ms. Long-dash - dot 
line describes M — R cq relation for neutron stars based on 
the FPS EOS, rotating at the maximum observed pulsar 
frequency. Curves for rotating neutron star models termi- 
nate at the low-mass end at the Keplerian (mass shedding) 
configuration, denoted by an open circle. 

structure, and can be described using slow rotation ap- 
proximation (Hartle 1967). Leading effects of rotation on 
neutron star structure are then quadratic in v, and are 
small as long as gravitational pull, acting on a unit mass 
element at the equator ~ G M / R 2 is significantly larger 
than the centrifugal force ~ Am 2 Rv 2 . However, with de- 
creasing mass of a neutron star rotating at fixed v, the 
rotational effects become more and more important. At 
the same M, the equatorial radius of a rotationally flat- 
tened neutron star becomes significantly larger than that 
of a static configuration. Importance of rotation increases 
with increasing equatorial radius and decreasing mass, and 
become decisive for v 2 ~ GM/(47r 2 i? 3 ). 

We performed exact calculations of stationary configu- 
rations of cold neutron stars, rotating uniformly at 100 Hz 
and 641 Hz. Numerical calculations were carried out us- 
ing exact 2-D equations for uniformly rotating configura- 
tions in general relativity. Equations of hydrostatic equi- 
librium were solved using numerical code, based on the 
pseudospectral method, and developed at DARC Meudon 
(Bonazzola et al. 1998, Gourgoulhon et al. 1999). Low- 
mass constant-^ sequences in the M — R plane are dis- 
played in Fig. 5. 

Effects of rotation are larger for the stiffer SLy EOS, 
which is most natural. At 641 Hz and M < 0.8 M Q , ef- 
fects of rotation are very large. The lowest-mass config- 
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uration corresponds to the mass-shedding limit. We get 
M^(641 Hz) = 0.61 Mq. For the FPS EOS, the effect is 
smaller, because it is softer than the SLy one, and there- 
fore low-mass neutron stars are more compact. We obtain 
M^f n s (641 Hz) = 0.54 M . Even at v = 100 Hz, rotation 
has a sizable effect, with minimum mass, reached at the 
mass-shedding limit, M^(100 Hz) = 0.13 M Q , larger, by 
nearly 40%, than that for static neutron stars. 

For the SLy minimum mass configuration rotating at 
v = 641 Hz the equatorial and polar radius of the liq- 
uid core are 11 km and 9 km, respectively. However, the 
crust is much more deformed than the whole star, hav- 
ing equatorial thickness of 6 km, four times larger than 
the polar one. This results in quite a large oblateness, 
-Rpoio/ficq = 0.6. The mass is concentrated in a liquid 
core, which contains 93% of neutron star mass. 

Minimum mass configuration, calculated for the SLy 
EOS at v = 100 Hz, has the (nearly spherical) liquid core 
constituting half of the stellar mass. However, the equa- 
torial radius of the liquid core is only about 6 km, to 
be compared with the 28 km equatorial thickness of the 
crust. Consequently, the oblateness of the minimum mass 
configuration is large, R po i e /R eci ~ 0.7. 

Our calculations show, that the rotating configurations 
close to minimum mass cannot be approximated by a ho- 
mogeneous rotating ellipsoid, like the model of Colpi ct 
al. (1991). The two quantities relevant for the Keplerian 
(mass shedding) limit, mass and equatorial radius, are de- 
termined by different parts of a low-mass neutron star. 
Namely, the mass is supplied by the liquid core, while the 
equatorial radius is defined by the outer edge of the ex- 
tended solid crust. Therefore, the approximation in which 
a neutron star is represented by a homogeneous (constant 
density - incompressible) spheroid (Colpi et al. 1991) is 
inadequate. 

It should be mentioned, that even at the mass- 
shedding limit at 641 Hz, rotating configurations are 
not susceptible to the viscosity-driven tri-axial instabil- 
ity which would imply emission of gravitational radiation. 
This secular instability sets in at kinetic to gravitational 
potential energy ratio T/W = 0.14. We get T/W = 0.061 
for the SLy EOS, and even less, T/W = 0.056, for the 
softer FPS EOS. Our conclusion is opposite to that ob- 
tained by Colpi et al. (1991), which results from the inad- 
equacy of the model used by these authors. 

One may ask, whether the values T/W = 0.05 — 0.06 
are above the threshold for the gravitational-radiation 
driven Chandrasckhar-Fricdman-Schutz (CFS) instability. 
The threshold values T/W ~ 0.04 - 0.05 for the m = 4, 5 
CFS instabilities, calculated for realistic EOS by Morsink 
et al. (1999), correspond toM^ 0.8 M Q . Extrapolation 
of these results to strongly deformed configurations with 
M ~ 0.6 Mq rotating at 641 Hz (Fig. 5) would be risky. 
Therefore, the problem of the CFS instabilities in rotat- 
ing low-mass neutron stars has to be clarified in the future 
studies. 

One has to point out basic difference between low-mass 
stationary uniformly rotating configurations and the static 



ones. Let us denote central density of the minimum mass 
configuration by /o c ,M min - In the static case, equilibrium 
configurations with p c < p c ,M min are secularly unstable 
with respect to radial perturbations. However, for v <; 100 
Hz, stationary configurations with p c < p Ci M min just do not 
exist. 

7. Summary and conclusion 

In the present paper we determined the properties of the 
minimum mass configuration of neutron stars, built of cold 
catalyzed matter. We used two recent EOS of cold cat- 
alyzed matter, which describe in a unified way both the 
crust and the liquid core, as well as their interface. For 
static neutron stars, we get a minimum mass of 0.09 M Q , 
coinciding with the value of M m ; n obtained three decades 
ago in the BPS paper. The value of M m ; n for neutron stars 
built of cold catalyzed matter is quite independent of the 
specific EOS used. However, the structure of the minu- 
mum mass configuration is sensitive to the EOS at p <, po- 
For the FPS and SLy EOS we get R(M min ) ~ 220 km 
and 270 km, respectively. Our minimum mass configura- 
tions contain a tiny liquid core with 2 — 3% of the stel- 
lar mass. The structure of the M m ; n configurations for 
the EOS considered in the present paper are substantially 
different from that obtained in the classical BPS paper. 
While BPS got a similar value of M m - m , the radius of their 
minimum mass configuration was only 160 km, and it was 
completely solid. These differences can be easily under- 
stood in terms of the differences in the EOS at subnuclear 
density. Namely, for p <J po, the BPS EOS is softer, and 
its crust-core transition takes place at significantly higher 
density than for the SLy and FPS ones. 

Within mass interval M min = 0.09 < M < 0.17 Mq, 
neutron stars are unbound with respect to dispersed 
56 Fe. In the similar mass range, neutron stars are less 
bound than carbon-oxygen white dwarfs of the same mass. 
However, all neutron stars are bound with respect to dis- 
persed hydrogen. 

Effects of rotation have been shown to be important 
for M m ; n , and for the low-mass neutron stars in general. 
For neutron stars rotating rigidly at P = 1.56 ms, mini- 
mum mass is determined by the mass-shedding limit: it is 
0.5 M and 0.6 Mq, for the FPS and SLy EOS, respec- 
tively. The minimum mass configuration at P = 1.56 ms 
is quite oblate, due to strong deformation of the crust by 
the centrifugal force. 
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